Diffusion of the electromagnetic energy due to the backscattering off Schwarzschild 

geometry. 
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Electromagnetic waves propagate in the Schwarzschild spacetime like in a nonuniform medium 
with a varying refraction index. A fraction of the radiation scatters off the curvature of the geometry. 
The energy of the backscattered part of an initially outgoing pulse of electromagnetic radiation can 
be estimated, in the case of dipole radiation, by a compact formula depending on the initial energy, 
the Schwarzschild radius and the pulse location. The magnitude of the backscattered energy depends 
on the frequency spectrum of the initial configuration. This effect becomes negligible in the short 
wave limit, but it can be significant in the long wave regime. Similar results hold for the massless 
scalar fields and are expected to hold also for weak gravitational waves. 



I. INTRODUCTION. 



Backscattering is a phenomenon that prevents waves from being transmitted exclusively along null cones. That 
aspect of waves propagation has been investigated for a long time for various wave equations (see, for instance, 
', @)- It has been established that solutions of the Klein-Gordon with nonuniform coefficients generically do exhibit 
■ backscattering jlj. This topic has been investigated in general relativity since early 1960's ( 0, [Q, i|, Q); a 
comprehensive bibliography can be found in @. The propagation of electromagnetic waves and of the resulting tails 
has been studied in the early seventies ( Q and j?J]) and recently by Ching et al. || in the context of the Schwarzschild 
spacetime and by Hod |)| in the context of the Kerr spacetime. The backscattering effect can be understood as the 
result of waves propagation in a nonuniform medium with a varying refraction index |Io| ] . 

In [ pd| has been assessed a classical aspect of the phenomenon that was not previously studied - the energy diffusion 
through null cones - in the example of a spherically symmetric massless scalar field propagating in the Schwarzschild 
geometry. The novel aspect of that work was a compact estimate on the magnitude of the backscattered energy in 
terms of the energy of initial data. 

This paper is dedicated to the investigation of propagation of electromagnetic fields in a background Schwarzschild 
spacetime. Similar to |]llf the main attention is focused on obtaining bounds on the backscattered fraction of the 
radiation energy, in terms of initial data. ^From the notional point of view the present paper parallels pd[ | , with three 
notable exceptions. First, the crucial technical points of the former work could have been applied only to spherically 
symmetric fields. In order to overcome this difficulty, the electromagnetic fields have to be split, with the extraction of 
; I ■ a known part which defines initial data. Then the standard expansion in terms of vector spherical harmonics leads a 
problem that can be tackled with methods applied earlier in pd| ]. Second, an energy inequality is proven. Third, this 
paper shows that the energy diffusion depends on the frequency of the radiation. An example of a dipole radiation 
allows one to characterize this quantitatively. The magnitude of the backscattering can be characterized as the ratio 
of the backscattered energy versus the initial energy of outgoing waves. This is vanishingly small in the short wave 
regime but it can be quite significant in the long part of the radiation spectrum. This kind of dependence on the 
frequency can be expected to hold also for higher multipoles. The scale is essentially set by the gravitational radius 
of the gravity source. All results of this paper hold true for any material sources of the Schwarzschild geometry - 
including stars, white dwarves, neutron stars and black holes, although the effects can really matter only in the two 
latter classes of objects. 

The order of the remaining parts of this paper is following. The next section defines notation, basic equations and a 
decomposition of the electromagnetic potential. The subsequent sections of this work deal only with dipole radiation. 
In section III is derived an energy estimate. Section IV is dedicated to the derivation of a bound, depending on the 
initial energy, of the backscattered part of the potential. Section V is devoted to the derivation of useful estimates of 
a pair of null-line integrals. In Section VI the equations are formulated in the language of characteristics. Previously 
found restrictions on the backscattered part of the potential allow one to estimate radiation intensities. Section VII 
brings an improved estimate of the backscattered potential, again basing on the method of characteristics. The next 
Section proves the main results - a bound on that fraction of the energy that can diffuse due to the backscatter off 
the Schwarzschild geometry curvature. Section IX shows that in the case of short-wave radiation the dipole radiation 
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backscatter is negligible. In contrast, in the long-wave regime the effect can be significant. Section X discusses how 
the effect depends on a distance and evaluates the exactness of the obtained criteria. The last Section presents a short 
summary and conclusions. 



II. FORMALISM 



ds 2 = -(1 - -j^)dt 2 + - \ m dR 2 + R 2 dti 2 , (1) 



Spherically symmetric geometry outside matter is given by a Schwarzschildean geometry line element 

2m N ,. 9 1 

where t is a time coordinate, R is a radial coordinate that coincides with the areal radius and dQ 2 = d0 2 + sin 2 ' 
is the line element on the unit sphere, < 4> < 2ir and < 8 < ir. 

As it concerns the electromagnetic fields, it is convenient to assume that the scalar component of the electromagnetic 
potential vanishes while the vector potential satisfies the Coulomb gauge condition. Using a multipolc expansion of 
the electromagnetic vector potential in terms of vector spherical harmonics one obtains [lj] 

+ = (l-^B + V*,. (2) 

\l/'s should be essentially two-index functions, *&im> (where M is the projection of the angular momentum), but since 
the evolution equation is <p— independent, the index M is suppressed. The variable r* = R + 2m\n(^ — 1) is the 
Regge- Wheeler tortoise coordinate. The backreaction exerted by the electromagnetic field onto the metric has been 
neglected in the present analysis. That is readily justified for any gravitational sources other than black holes. In the 
case of a black hole this approximation holds true some distance away from its horizon fl3|| . 
Consider a set of functions of the form 

s=0 

where the functions \I/j s are given by the recurrence relations 



1(1 + 1) 
df*li = — L -2- L *lo 



(s(s + l)-l(l+l))V ls -2m(s 2 -!)*,(,_!) . (4) 



In is shown a dipole solution of this type O]. In the Minkowski space-time (m=0) ^1 solves (||); it represents a 
purely outgoing electromagnetic radiation. 

Let a function 'J/ be given by (||) and ([|) and assume that (for space-like sections with t > Q) its support is compact 
and located entirely in the vacuum region outside some radius a > 2m, i. e., outside the Schwarzschild radius. Let 
initial data of a solution of (||) concide with ^>i at t = 0. Thus initially ^1 is a purely outgoing partial wave. It 
should be noted that the assumption that initial data are (initially) purely outgoing is made in this paper only for 
the sake of clear presentation. The propagation of electromagnetic waves is a linear process as far as the backreaction 
can be neglected. Therefore the propagation of the initially outgoing radiation (or even of a fraction of the outgoing 
radiation) is independent of whether or not the ingoing radiation is present. 

It will be convenient to decompose the sought solution ^/(r*,<) into the known part and an unknown function 

Si 

^1 = ^1 + Si. (5) 

Initially Si — doSi — 0. A similar splitting is done in Q, who then seek a series expansion of Si. This will be avoided 
in the paper, in favour of finding a number of estimates of Si that would provide the needed information about the 
backscattered part of the radiation. 

In the rest of this paper will be considered only dipole radiation, ^1. Consequently, all angular momentum related 
subscripts will be omitted. 
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III. AN ENERGY ESTIMATE 



The dipole term constitutes the most important part of the electromagnetic radiation. Assume dipole-type initial 
data 



*(x(R)) = -d r *f(x(R)) + 



R(r*) ' 



(6) 



with the initial support (a, oo)) of a C 2 -differentiable / and x(R) = r* (R) — r* (a) . The differentiability of / guarantees 
that the initial energy density is continuous and vanishes on the boundary a. 
Lemma 1. Define I a e (R) 



Ia,e(R) 



dr 



f 2 (x(r)) 

r i+2e ■ 



and 



(7) 



(8) 



where < 2e < 1. Then for a > 2m(l + 1/y/l + 2e) the following inequality holds 



f3a(R) 



(*)' 



(9) 



ea 2e (l + 2e)(l-^f) 2 -^' 
Remark. The integral (3 a {R) is bounded above by the electromagnetic energy En(t)/(4:Tr) defined later. Therefore 



Ia.e(R) < 



E a (R,t) 



2 c. 



Airea 2e (1 + 2e)(l - ^f) 2 - ^ ' 



Proof. 

Notice that 



P(x(r)) 



dsd. 



f(s) 



ds 



i 2 m 



cis- 



W 2m/ 
1 

s s 3 

2 

+ 2 



< 



ds- 



2mf 



>(1 -if) 



(10) 



the second inequality follows from (A + B) 2 < 2A 2 + IB 2 . The factor 1/(1 - 2m/ s) that appears in the integrands 
can be bounded above by l/t] a , where 



2m 
a 



Va = 1 

Subsequently, the use of the Schwarz inequality and simple integrations yield 

P(x{r)) < 

1 



(11) 



2/3 a (r){r-a) , 8m 2 / a , £ (r) / 1 



r? 2 (l-2e) Va 1 



l-2e „l-2e 



(12) 



The insertion of (fL2|) into the integral of ua) gives 
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dr- 



1 



.2+2e 



2/3 a (r)(r-a) , 8m 2 /„, f (r) 



'7a 



r?2(l-2e) Va 



2 c 



,,l-2r 



(13) 



/3 a (r) and I a ,e{ r ) are nondecreasing functions, therefore taking them in front of the appropriate integrals would not 
make the corresponding terms smaller. Straightforward integration of the obtained expressions yields 



Va 

2J a , 6 (fl) / 2m ^ 
1 — 2e V a — 2m./ 



2c 



1 + G 



l+2e 



+ 



1 



l + 2e 



:i-(^ 



l+2e 



RJ 



1 , /a 



(14) 



One should note that the expression inside the bracket in the first line can be estimated as follows 



1 / / a \ 2< \ 1 / /a 



l+2c 



< 



1 



(1-f- 

e(l + 2e) V Vi? 



while the expression inside the bracket in the second line is bounded above by 



1 - 2e 
2(1 + 2e) 



l+2e 



Eq. ( |l4| ) can be now written as 



IaAR) < 



( \ 2e 

(3 a (R)(l - ) . 2m \ 2 / a , £ (r) 

r? 2 a 2e e(l + 2e) + Va-2m/ 1 + 2e ' 



Rearranging Eq. ( |15[ ) so that the two terms with I a ^(R) are on the left hand side, one obtains 



2m 



1 



a -2m J 1 + 2e 



< 



/3 a (fl)(l-(^) 2£ ) 
77 2 a 2e e(l + 2e) 1 



(15) 



(16) 



this gives the postulated bound of I a ^(R) if a > 2m(l + 1/^/1 + 2e), as assumed above. 

The obtained formula is not exact, but with the appropriate choice of / and e the error is small. Take for instance 
f = C within (a + ai, b— cii,bi << a, 6 >> a (which obviously means that b — a >> a\, &i), and let / be smoothly 
joined to zero outside (a, 6) by some intermediary functions. Under those conditions, a direct calculation gives 



IaA b ) „ 3 

A, (6) ~ (3 + 2e)a 2e ' 



(17) 



as compared with e ( 1+ 1 e - )a2g that follows from (16). If e « 1/2, then the exact result differs by less than 25 percent 
from the bound of ( |l6| ) . Later on will be used the value e = 1/8 (which appears to be more economical in subsequent 
calculations), in which case the above estimate deteriorates significantly. The exact value of ^''(ff is then roughly 15 



% of that predicted by (fffl). 



IV. ESTIMATING S. 



S is initially zero and its evolution is governed by the following equation 

6m/ 



(-d 2 + d r 2 ,)<5 = (l-^) 



R 4 



(18) 



Define t ) - a null geodesic that originates at (R, t) and is directed outward. If a point lies on the initial hypersurface, 
then f (fl)0 ) = f R . By f ( Ro>to ),(R,t) will be understood a segment of f (iWo) ending at (R,t). 
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Later will be needed a following bound. 

Theorem 2. Let the support of initial data be (a, 6), b < oo and let TrqIrj) be the outgoing null geodesic from 
(Ro,t = 0) to (R,t). Then 



7? 



where 



Ci = 



6^2 



^ /2 (l-e)Ve[(l + 2^-^] 



(19) 



(20) 



and 7/^ = 1 — 2m/ R. 

Proof. Define an energy H (R, t) of the field 6 contained in the exterior of a sphere of a radius R, 



H(R 



(21) 



One can easily show that 

(d t + d;)H{R,t) 



2m , / 9n<5 



(1-— )(- 



12m / dr<9 <5^j < 



12m / drdoS- 



(22) 

the inequality follows from omission of the nonpositive boundary term. The right hand side can be bounded further 
by 



12m 



dr(d Sf 



1/2 



00 / 2 ]l/2 
R r J 



due to the Schwarz inequality. That in turn can be bounded by 



12m 



H 1 ' 2 



dr 



r 4+2e 



1/2 



(23) 



(24) 



The second integral in (|2J) can not increase along outgoing null directions, therefore is bounded by initial values, 
dr^lvzt < Jr dr^i^ = lR a ,e( R )- Since I R a ,e{oo) < J , 6 (oo), one arrives at 



(d t + d;)H(R,t) 1 ' 2 < ey^M- 



^TaR 2 - 



(25) 



The integration of ( |25| ) along Tr q ^r^ yields, replacing J a c (oo) by its bound expressed in (g), 

1 



y/H(R,t)< 



6m y//3 a (c 



1 



1 



( ol-c E>l-e ) 



« £ 11 l/ 2 {l-e)\je{{l + 2e) V l~^Rl^ i? 1 " 
mC\ , 1 1 , 



(26) 



Notice that initially 5 vanishes and that its propagation is ruled by a hyperbolic equation. Thence at any finite time 
t the support of S is bounded. Therefore 



\5(R)\ 



R 

R 



< 



1 n1/2 



<^2\i/2 



( i*) 1 ( (^-^) 2 ) 



< 



1 

\AR-2m 



{2H) 1 ' 2 {R) 



(27) 



5 



Inequalities ( |2q ) and (27) yield the bound of Theorem 2 in the case when b — oo. 

Let the initial data be of finite support (a,b). Define a region fi& consisting of points (R > b, t) acausal to (b, t = 0). 
The energy H(b(t),i) obviously vanishes for any point (b(t), t) located inside fib- In this case the inequality of Theorem 
2 can be stated as follows: 



— < mCWm (^7 - J^h (28) 

In what follows it will be assumed that the initial data have compact support located in an annular region (a, b). 



V. ESTIMATES OF TWO (NULL) LINE INTEGRALS 

In analogy with fWt) defined earlier, let L( fljt ) be a null ingoing geodesic that originates at (R,t). Tm t —o) will be 
shortened to T R . A segment of Tr Rljtl \ connecting (Ri,ti) with (R, t) [t\ <t,R\> R) will be denoted as Tr Rlttl \r Rjt \. 

Let a point (R, t) be an intersection of an ingoing null geodesic F Rl with an outgoing null geodesic r a . Let (r, r), 
r > R be a point of T Rl ^ R t ^ and define (Ro(r), t — 0) as a point of the initial hypersurface such that r^nT^ = (r, r). 
Fixing a and Ri, one can view i?o as a function of r; obviously Ro(R) = a while Rq(Ri) = R\. On the other hand, 
fixing only a and viewing R\ as a function of R one has Ri(a) = a; this will be used in the forthcoming proof. 

One can prove 

Lemma 3. Under above conditions and if R\ < b, the line integral along a null segment geodesic F^^^ is 
bounded above, 



Ri 



dr- 



Rq 



1 /, R 
< - In — 



ln( 



b — 2m , 
2m 



a 



(29) 



Lemma 4. Under the above condition but with the initial point (Ri,s) (s > 0) of the null geodesic segment 
^(Rx,s),(R,t) lying on f b (Fig. 2), one can prove 



Hi 

dr — 
RorJl 



Ri 



1, ,6- 2m, 
= < - In 5- 

2m 2 a — 2m 



(30) 



Proof of Lemma 3. 

Let r be a radial coordinate of a point lying on the intersection of T Ro and T Rl (Fig. 1). One finds that the areal 
distances of three points (Ro(r),0), (r,r) and (Ri,Q) satisfy following 



Ro(r) = 2r- Rx + 2m In 



(r - 2m) 2 



(i?o(r) -2m) (Hi -2m) 



(31) 



That implies 



dR = 2 



2 zrn 



i 



2m 



dr. 



(32) 



Replacing r by Rq in the integral of (E9J) one obtains 



Hi 



dr- 



r-Rn 



RorJl-^ 



Hi 



dRo 



1 



Ro - 2m 



1 ^-2rn * 

2 v a - 2m ; i? 



Hi 



Hi 



In- 



I in 1 



2m 
Hi 



dr I \ 1 

_ 2m Vi? r 

r 



< 



(33) 



G 



Next, one can show that R\ > 2R — a. Indeed, assuming that a is fixed, one has from ( |32| ) > 2; since the initial 
condition is Ri(a) = a, the conclusion follows. 
Taking into account R\ > 2R — a one gets 

R . R / R . , 

In ■ <ln — — <ln\— (34) 

VoTTT y/a{2R- a) V a 

Replacing Ri by b in the last term of (|33| ) and inserting ([m|), one arrives at the first of conjectured inequalities, 

In order to prove Lemma 4 one should start from relation between areal distances of four points (r, r) , (i?o (r) , 0) , 
(R,t) and (o,0) (see Fig. 1): 

2(r-R + 2mhx^^) = Rq- a + 2mln R °~? m ; (35) 
V R— 2m ) a — 2m 

The variable r ranges from R\ > b to i?. Fixing a and R, one again obtains 

^ 2m 

<*flo = ^Y^dr. (36) 

r 

A straightforward calulation, in which dr is replaced by dRo, shows that 
Since i?i > R, one immediately obtains (fj[)|). 



VI. AN ESTIMATE OF THE AMPLITUDES BACKSCATTERED INWARD 



Define the intensity of the backscattered radiation that is directed inward 

1 



Eq. @ reads now 



h-(R,t) 



2m, 



1 



-2m (^o +d r *)5. 

R 



(-d + d r .)[(l-—)h- 



2m, 



6mf 



(38) 



(39) 



The integral form of (B9J) reads, 



(1-— )M*>*)= / 



/ dr 


2 ^ ^ 6mf 


IRx 





h-(Ri, s); 



(40) 



here the integration contour coincides with a null ingoing geodesies ^(R 1 . s ),(R,t)- s ) nes on the initial hypersurface 
(s = 0) if i?i < b; thus h-(Ri,s = 0) = 0, since the initial data are entirely outgoing. If i?i > b then (i?i, s) € IV( fll . s ); 
also in this case h-(R\, s) = 0, because ^b,Ri constitutes the outer boundary of the outgoing impulse. In either case 
the radiation amplitude satisfies the integral equation 



(1 



2m 



f R 

h-(R,t)= / 
Jr.! 



dr 



2 s , 6m/ 
-r<5 H r- 



(41) 



The second term is bounded above by 
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L 



Ri tl 



6m 



1/2 



cir- 



ri 



.4+26 



1/2 



dr- 



4-2e 



6m 



V3 - 1 



Hi j2 



1/2 



4+2e 



r 



#( 3 / 2 )- 



3-2e 



< 



6m 



/ 2 



1/2 



1 



r 4+2e I R(3/2)-( 



1 - 



,3-2e 



5 3-2 £ 



(42) 



where the second line follows from the Schwartz inequality and the last inequality is due to the fact that R/Ri > a/b 
(Appendix A). 

In order to find the integral from the last line of of fl42| ) it is useful to project it onto the initial data surface, along 

l- 



R r 



^4+2e 



clR Q - 



2 m 



outgoing null geodesies T Ro ^ r ^ T y Notice that dR = 2 ^ R §^ dr - see (|2|). The J 2 / r 4 + 2e term cannot decrease during 
this projection. One arrives at 



Inserting the energy estimate of Lemma 1 into ( ff3| ) one gets finally 

mC 2 



(R ) < IaARl) 



2r] a 



1 -<?)-. 



Here the constant C2 is given by 



3^2(1 



^ /2 V /e[(l + 26) ?? 2 -^](3-2 e ) 



2 _ 4m 2 1 



(43) 



(44) 



(45) 



The ^-related term of ( [4l| ) is bounded, due to ([lj]), by 



2mCi 



(46) 



Here r > Rq and r, i?o G f fl 0i ( r , T )- Thus l/(r £ i?Q e ) < 1/Rq. Therefore the expression of ( fig ) is bounded above by 



2mCi 



VfUb)- 



yr- TIF f Rl A 1 / 1 1 



i? 3 / 2 - e 

Results of Lemma 3 and 4 lead now to a pair of estimates. If R\ < b then 
and if i?i > 6 (in which case (i2, i) € IVa.^)) then 



2 / 1 dr^ <mC lv ^(6) 



ir4 < mdV^) 



e R 3/2 



Va 



(47) 



(48) 



(49) 



In summary, the radiation amplitude is bounded above by 
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(l-—)\h_(R,t)\< 



C 3 

e R 3/2- 



C 4 + d In 



66(i?i(i?) - b) + RQ{-R 1 {R) + 6) 



(50) 
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where 0(6 - R±) = if 6 - Ri < and 9(6 - Ri) = 1 if b — Ri > 0. The constants C 3 and C4 are defined by 

C 3 = my/^JlTfy*, 



C i = C 2 + C 1 lD.—. 



(51) 



VII. REFINING THE BOUND ON 5 



Equation (Bq) can be written in the integral form 



S(R,t) 



drh- + 5{R ), 



(52) 



where the integration contour coincides with F(r q ) jr ,t) an d i?o is a point of the initial Cauchy slice defined earlier. 

Since initially 5 vanishes, one has 5(R,t) = drh-. It becomes clear in Sec. VIII that one needs to bound S(R,t) 

only along r a; oo; in what follows is always meant this situation. Define R(b) = (R : Ri(R) — b) (see Fig. 2). Inserting 
the bound of (|50| ) (but notice that (|o"o| ) bounds 7] r \h{r)\, not \h(r)\ itself), one obtains 



\W)\ < 



r R dr f H In- 



c 3 



dQ(R- R(b))( I dr 



3/2-6 



In 



dr 



(53) 



The integrand of the second integral is nonnegative, therefore extending the integration up to R(b) can give only a 
bigger quantity. Thus one gets, after elementary integration, 



\S(R)\ < 
2C 3 



Va (l-2e)a^ 
1/2-6 



2Ci 
1 -2e 



1 - 



i?(6) 



1/2-1 



i?(6) 



In 

Dropping out the negative term — 9(i? — R(b)) In - 1 — g- ■ 
-R(6)) In - < In -ht^t, one arrives at 



e<«-«<»„ h A (l _(M r 



(54) 



6 /-R(b)\ 



1/2-e 



and taking into account that — ln-^^ + Q(R 



2C 3 



Va (l - 2c)oVa 
2Ci 



C 4 1 



R 



1 - 2e 



1 - 



1/2-1 



Define 



One can show (see Appendix B) that 



R{b)J J \R(b) 
k = (6 — a)/a. 



a + b s ntu\ s a + b 

— run < R(b) < — — ; 



R(b) 



(55) 



(56) 



(57) 
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the equality is achieved in the Minkowski space-time (m=0). Since b = a + an, one has R{b) > a + r) a n/2 or, defining 

a=^, (58) 



R(b) > a + an. The insertion of the above into (|55| ) yields 

C 3 



where 



C5 = C\ - 



77 a (l - 2e)aV2 



«('-(!) 



1/2— e 



C 5 



In 



2Ci 



^1 + a/t^ 



1 - 



1/2-6 i_ 2e V (I + k/2) 1 / 



(59) 



(60) 



This estimate gives a better control over the asymptotic behaviour of 8 than the former one, (|27|), by a factor 
In particular, now S 2 / R 2 is known to be integrable. This integrability will be exploited in the next section. 



VIII. BOUNDING THE RADIATION ENERGY LOSS 



The energy En(t) of the electromagnetic field W contained in the exterior of a sphere of a radius R reads 



E R (t) 



2k I ilr\ 



(W 2 

-y 2 m 



(1 



2m, 



2(*) : 



(61) 



Let the initial data be as specified hitherto, ^>(t = 0) = \& and c^o^^ = 0) = c'o'I' for some \&; thus they vanish 
outside an annular region (a, b). Define = E a (0) as the energy of the initial pulse. If initial configuration is purely 
outgoing, then <9 t * = -9 r »# + fd r *(l/R). 

Let an outgoing null cone C a originate from (a, 0). In the Minkowski spacetime the outgoing radiation contained 
outside C a does not leak inward and its energy remains constant. In a curved spacetime, however, some energy will 
be lost from the main stream due to the diffusion of the radiation h- through C a ■ Most of the backscattered radiation 
will fall onto the center of the gravitational attraction. The forthcoming theorem gives a bound on the amount of 
diffused energy. 

Theorem 3. Under the above assumptions, the fraction of the diffused energy 8E a jE h a satisfies the inequality 



SE a /2m\ 2 l ~ (T+^y 



El 



\ a / 



C\ (Tja 



2e 2 



C? 



(l-e)Vl6 (1 - 2e) 2 (3- 2e) 
where G\ - C5 have been defined earlier and 

VaC 2 / 



(l-2e) 2 



C fi 



2C4C5 



(1 - 2e)(3 - 2e) 



Co 



G 



2 l+K 
1 + QK 



1 



+ 



16(1 -e) V(l + a K ) 2 " 2e (1-e) 2 

Ciln(l+K/2) ^ 1ti 1+k x _ 



CiC 4 ?7 a /^+ c 4 (i-e) U n T+S^ 



16(1 -e) 



(1 



\2-2c 



(1+k/2)*~^ 



Proof. The rate of the energy change along C a is given by 

{d + d r .)E a = 



-2?r(l 



2m N 
IP 



(1 



2m 



2m, 

-Mi -^5 



0«* 

2 



i? 2 



* 2 



2mJ f \ 2 /- 



J? 



J? 2 



i? 2 V 



(62) 



(63) 



(64) 
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The functions / and W are assumed to vanish on the null cone C a . Therefore VP = 8, <9r 1 J / = d^S and dt^ = dtS on 
C a . In such a case the rate of the energy change reads 



(d + d r ,)E a = 

2m, r 
~R 



-2tt(1 



((1 



2m x 



2^2 



(65) 



The energy loss is equal to a line integral along T a , 



fOO 



2n 



dr 



- — 

2m,, , 25 2 
(! - —) h - + — 



(66) 



The derivation of ( p2\ j requires the use of estimates ( |50| ) and (59). The calculation of the 5- related part of the right 
hand side of fl66|) is straightforward and it yields 



4tt 



(v n 
— 
a 



2m\2 1 



r2 2e 2 
4 (3-2e)(l- £ ) 



C K + 2C4C5 



l-2e 
3- 2e 



(67) 



In order to bound the contribution coming from the backscattered radiation amplitude h— one needs the estimate 
(fiol) . A straightforward calculation shows that 



/■°° 2m 
2n / d r (i - fH)h 2 _ < 

J a ' 



2m y 1 ~ (i+k)^ 
a J 2^(1 -e) 



1 



Cf 



2 y 



hr 



J?(6) 2(1 -e) 



i?(6) 2(1 - e) 



In- 



(1-W)J + 
(1-2 In " 



i?^(6) y 4(1 -ef 



(1-2/) 



(68) 



where y = (a/ R(b)) 2 ~ 2e . Neglecting the negative term with In... and using the bounds of Appendix B on b/R(b), one 
arrives at a bound that in conjuction with ( |S7| ) proves Theorem 3. 

Remark. The above estimate depends on the parameter e, that should be chosen in such a way as to optimize the 
bound. The exact value of the optimal e depends on a and k, but the value e = 1/8 is proven to yield satisfactory 
estimates. 



IX. DEPENDENCE OF BACKSCATTER ON THE FREQUENCY OF WAVES 

The coefficients C4 - Cg appearing in Theorem 3 change with k but remain finite in the whole (0, 00) range of 
possible values of k = (b — a) /a. In the case when the support of the initial radiation is very narrow, i. e., k « 1, 



1-- 



then the coefficient — ( 1 +"T^ r ^ K j sucn a case one obtains that 

1 * c '(v) 2 - «•» 

where C is a constant. In the limit k — ► the ratio ^ L becomes 0; the backscattering is negligible in the case of initial 
pulses of electromagnetic energy that are very narrow. And conversely, the bound becomes bigger with increase of 
the width of the radiation pulse. The physical meaning of that can be deduced with the help of the Fourier transform 
theory. The similarity theorem ( (T5) ) states that compression of the support of a function corresponds to expansion 
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of the frequency scale. If a support of initial data is made narrow, then the wavelengths scale of the pulse extends 
in the direction of short lengths. Therefore the message behind the obtained results must be that high frequency 
radiation is essentially unhindered by the effect of backscattering and that long waves can be backscattered. 

This dependence of the backscattering on the wave length has been in fact observed in the numerical investigation 
of the propagation of pulses of scalar massless fields . In this case halving of the length has led to a similar decrease 
of the fraction of the diffused energy. 

In the case of a black hole or a neutron star the scale is set essentially by the Schwarzschild radius Rs — 2m; waves 
with lengths much shorter than Rs are not backscattered, while waves of lengths ~ Rs can reveal quite a strong 
effect. Moreover, one can show that the (2m/ R) 2 dependence of the effect implies that most of the energy diffusion 
occurs in regions that are not very far (as compared to the Schwarzschild radius) from the center. 

In order to exemplify the above remarks, consider the diffusion effect in following two cases. Assume the same 
location a = 4i?g, of both radiative dipoles and 

i) k = 1/8 (i. c., the fundamental wavelength Rs) for a pulse I; 

ii) k = 1/128 (i. e., the fundamental wavelength Rs/8) for the pulse II. 

In the calculation that is reported below, e is chosen to be 1/8, in accordance with the remark ending the preceding 
section. Then in the case I one obtains 5E a / E h a < 0.37, while in the case II (of shorter waves) one gets 8E a jE h a < 0.004. 
The evolution equation (|39| ) can be written in another form as 



(d + 9 r .)(0 



2m 
~R 



(1 



2m x 
~~R' 



2 6m/ 
~B* + ~R r 



where 



h+ (R,t) = -^-^ (-do + d r , ) (6 + /) . 



(70) 



(71) 



is the intensity of the outgoing part of the radiation. The inequality ( |l2| ) can be written as follows, applying Lemma 
1 and the remark following it, 



(72) 




here C is some constant. The integration of (|70| ) along a null geodesic f a yields now 



2m, 



2m, 



(1 - —)h + (R(t), t)-(l- —)h + (R(0), t = 0) cx I 1 - 



2c 



(73) 



where the proportionality constant depends only on e, 2m /a and the initial energy E h a . Fixing the energy E^, one 
notices that in the regime (b — a)/a << 1 the right hand side of (|73| ) is essentially zero. Thus the product (1 - 
is constant. In this case one clearly sees the manifestation of the redshift - the rescaling of the amplitude h + 



h + (oo) = T] a h+(a,t = 0). 
See also a discussion of that fact in a massless scalar field theory jrl) . 



(74) 



X. DISTANCE DEPENDENCE OF ENERGY DIFFUSION AND SHARPNESS OF THE ESTIMATES. 

The bound of Theorem 3 depends on the source location - it contains, among other factors, a square of the factor 
2m/ a. Thus the bounds in question decrease with the increase of a. The dependence on the distance can actually 
be much stronger. In order to see this, consider the dipole radiation II of ii), described in the preceding section, 
but located at a = 4m (instead of a = 8m, as assumed formerly). One obtains that now 8E a jE h a « 0.77, instead of 
the former bound 0.001. Numerical results concerning the propagation of massless scalar fields also show that the 
backscattered energy decreases rapidly with the increasing distance jllj . 

It is of interest to establish how accurate the final estimate is. Most of the inequalities derived in this paper are 
sharp, in the sense, that one can find examples that saturate them. Thus, for instance, the two null-line integrals of 
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Section V are estimated sharply (the inequalities saturate in Minkowski space-time). Similarly results of Appendices 
A and B are also exact; again, the inequalities become equalities in Minkowski space-time. The energy estimate of 
Section III is not sharp; but the "loss of sharpness", to say, can be less than 25 % (see the final remark in Sec. 
III). The main source of unsharpness is the omission of negative terms in a bound on 8 (Sec. VII) and in bounds of 
diffused energy in Sec. VIII; but that becomes insignificant with the decrease of K, i. e., when the width of the pulse 
becomes small in comparison to the Schwarzschild radius. On the other hand, the combination of two exact steps can 
be associated with some loss in the accuracy. 

Taking this into account, it is quite likely, that in the case of sources characterized by k < 1 the bound in question 
gives an order of the diffused energy. On the other hand, if k » 1, then the bound of Theorem 3 becomes very 
inaccurate, with 



where C ~ 10 . As mentioned before, the main source of unsharpness is the omission of negative terms in a bound 
on 5 (Sec. VII) and in bounds of Sec. VIII. A more accurate treatment would significantly improve the estimates in 
the long wave regime. 



The main result of this paper, Theorem 3, states that the dipole energy diffusion due to the backscattering depends 
on the square of 2m/ a, where m is the mass of the gravitational source and a is a location of the pulse of radiation. 
Sections IX and X show that the high-frequency radiation essentially is not backscattered, but that the low-frequency 
radiation can manifest a significant diffusion effect. The last statement is best described in terms of the dimensionless 
parameter k = R$/X, where A is the fundamental radiation length. If k >> 1 then the backscattering is negligible, 
but if k fv 1 then it can be significant. The above results demonstrate that the effect becomes negligible at distances 
much bigger than the Schwarzschild radius of a central mass. That rules out most stars as objects that can induce 
observable backscattering effects. For a star of a solar type and A ~ Rg, for instance, the ratio ^g- can be at most 

10~ 20 . In the case of white dwarves and A ~ Rs the bound J75| ) gives ^ L < 10~ 8 . For long-wave radiation the 
bounds are bigger - the effect even looks as marginally relevant, for white dwarves, when -J^ 1 - ~ 10~ 2 . However a 
sharper estimate would lower that by several orders. 

On the other hand two astrophysical compact objects, neutron stars and (most likely) black holes, are not excluded 
as objects of interest. 

The backscattering would damp the total luminosity produced in accretion disks that exist in vicinities of compact 
objects, but since the most efficient regions of the disks are located at a distance of (at least) several Schwarzschild 
radii, the effect would be probably weak. More relevant can be "echoes" - aftermaths of violent flashy eruptions, 
produced by a part of radiation reflected from a close vicinity of a horizon of a black hole. Numerical calculations 
done in the massless scalar fields propagation suggest that the amplitude of the reflected radiation can constitute up 
to 20 % of the incident one, assuming that the length of the wave is comparable to the Schwarzschild radius of a black 
hole. 

The results of this section can be in principle generalized into the case of higher order multipoles. The key point 
would consist in showing analogues of the energy estimates of Sec. Ill that would bound the higher multipole moments. 
That should lead to a variant of Theorem 3 valid under reservations similar to those expressed earlier. 

An analysis similar to that of the present paper can be repeated also in the case of a weak gravitational radiation 
produced in disks rotating around Schwarzschildean black holes. The conclusions concerning the fraction of the 
diffused energy can be similar. 

Acknowledgements. This work has been suported in part by the KBN grant 2 P03B 010 16. The author is 
grateful to Niall O' Murchadha for reading of the manuscript, many discussions and valuable comments. Thanks are 
also due to Peter Aichelburg for a useful remark and to Krzysztof Roszkowski for his help in editing of this paper. 



[1] J. Hadamard Lectures on Cauchy's problem in linear partial differential equations, Yale University Press, Yale, New Haven 




(75) 



XI. DISCUSSION 



1923; 



13 



[2] C. Misner, K. Thorne, J. A. Wheeler, Gravitation, Freeman, San Francisco, 1973. 

[3] B. DeWitt and R. Brehme, Ann. Phys. 9, 220(1960); W. Kundt and E. T. Newman, J. Math. Phys. 9, 2193(1968); 
R. G. McLenaghan, Proc. Camb. Phil. Soc. 65, 139(1969); W. B. Bonnor an M. A. Rotenberg, Proc. R. Soc. A289, 
247 (1965); J. Bicak, Gen Rel. Grav 3, 331(1972); R. Price, Phys. Rev D5, 2419(1972); G. Schaefer, Astronomischer 
Nachrichten 311, 213(1990); R. Price, J. Pullin and Kundu, Phys. Rev. Lett. 70 1572(1993); C. Gundlach, R. H. Price 
and J. Pullin, Phys. Rev. D49, 883(1994); ibid, 890(1994); L. Blanchet, Gravitational Radiation from Relativistic Sources, 
in Relativistic Gravitation and Gravitational Radiation, eds J. A. Marck and J. P. Lasota, (Cambridge University Press 
1997); L. Blanchet, Class. Quantum Grav. 15, 113(1998); 

[4] J. M. Bardeen and W. H. Press, J. Math. Phys. 14, 7(1973); also p. 868 in C Misner et al. in |j. 

[5] W. B. Bonnor and M. Piper, Class. Quantum Grav., 15 955(1998). 

[6] L. Blanchet and G. Schaefer Class. Quantum Grav. 10, 2699(1993). 

[7] M. A. Rotenberg J. Phys. A4, 617(1971); W. E. Coach and W. H. Halliday J. Math. Phys. 12, 2170(1971); A. Papapetrou 
J. Phys. A8, 313(1975). 

[8] E. S. C. Ching, P. T. Leung, W. M. Sueng and K. Young, Phys. Rev. D52, 2118(1995); Phys. Rev. Lett. 74, 2414(1995). 
[9] S. Hod, Phys. Rev. 60(2000), 104053; Phys. Rev. Lett. 84, 10(2000). 
[10] see U for a relevant transformation of the equations. 

[11] E. Malec, N. O'Murchadha and T. Chmaj, Classical and Quantum Gravity, 15, 1653(1998). 
[12] J. Wheeler, Phys. Rev. 97, 511(1955). See also p. 868 in C. Misner et al. in §. 

[13] The condition for the smallness of backreaction is of the type shown in Jllj for the massless scalar field. 

[14] Incidentally, an expansion of Q must be divergent. Its covergence would essentially demand that |j"(ij(r))~ 1 j°° g{s)ds\ < 1, 

which cannot be true for functions g of compact support. For another argument see Bardeen and Press in H|. 
[15] R. Bracewell, The Fourier Transform and its applications, McGraw-Hill Book Company 1965. 
[16] P. Regucki, MSc. thesis, Institute of Physics, Jagellonian University 1999. 



14 




15 



XII. APPENDIX A 



Lemma A. Let (R,t),Ri £ T Rl j R<t \, R\ > R and R > 4m. Then 

R a 

K * 1 (76) 

Proof. 

There are two separate cases that need to be considered. 

i) If Ri lies on the initial hypersurface then R\ < b and R > a and the inequality follows immediately. 

ii) If (i?i,s) G ^b,(Ri,s)- I n this case one has 

2(R 1 -R + 2m\n Rl ~ 2m ) = b - a + 2mln^—^- (77) 
V R— 2m J a — 2m 

Define X = R\ — R. One finds from @ that 

provided that R < 4m. Thus X/ii is a non-increasing function which means that R/Ri is a non-decreasing function 
and R/Ri > R(b)/b. Since R(b) > a, one arrives at the postulated inequality. 

XIII. APPENDIX B 

Lemma B. Define k = (b — a) /a < 0. Define (R(b), t) as the the intersection point of Tj, and T a . Then 

a + b a + 6 , , 

^--m K <i?(6)<^-. (79) 

Proof. The relation (|3l| ) (see the main text), with R\ = b, r = R(b) and i?o = a, can be written as , 

a = 2R { b)-b + 2mlJ ( *f ~ ^ A (80) 
V (a — 2m) (6 — 2m) / 

We will treat ( |8C| ) as a relation between 6 and i?(&), with fixed a. Obviously R(b) = b = a when b = a. One easily 
finds that 

d b R{b) = i^iSL . (81) 
2 ?7b 

Notice that J?(6) < 6. Thus d b R(b) < 1/2. On the other hand R(b) > a. Thus d b R(b) > (l/2)rj R(b) > (l/2)rj a . The 
use of those two bounds on d b R{b) and the initial condition R(a) = a immediately imply the Lemma. 
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